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Models with dark energy decaying into dark matter have been proposed to solve the coincidence 
problem in cosmology. We study the effect of such coupling in the matter power spectrum. Due to 
the interaction, the growth of matter density perturbations during the radiation dominated regime 
\Q • is slower compared to non-interacting models with the same ratio of dark matter to dark energy 

' today. This effect introduces a damping on the power spectrum at small scales proportional to the 

strength of the interaction, c 2 , and similar to the effect generated by ultrarelativistic neutrinos. 
The interaction also shifts matter-radiation equality to larger scales. We compare the matter 
power spectrum of interacting quintessence models with the measurments of 2dFGRS. The data 
are insensitive to values of c 2 < 10 -3 but strongly constraints larger values. We particularize our 
study to models that during radiation domination have a constant dark matter to dark energy ratio. 
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I. INTRODUCTION 



Observations of high redshift supenovae , temperature anisotropics of the cosmic background radiation 0, , 
matter power spectrum 0, and the integrated Sachs- Wolfe signal @ indicate that the Universe is currently 
^sO • undergoing a phase of accelerated expansion Q- A cosmological constant, A, is frequently invoked as the most 
natural candidate to drive this acceleration. However, this choice is rather problematic. First, the observed A value 
falls by many orders of magnitude below the prediction of quantum field theories |]J . Second, it is hard to understand 
why precisely today the vacuum energy density is of the same order of magnitude than that of matter. This remarkable 
fact, known as the "coincidence problem" 0, lacks of a fully convincing theoretical explanation. 
' Models based on at least two matter components (baryonic and dark) and one dark energy component (with a 
high negative pressure) have been suggested to explain the accelerated rate of expansion and simultaneously alleviate 
the coincidence problem 0, [n|. Further, coupling between dark matter (DM) and dark energy (DE) has been 
suggested as a possible explanation for the coincidence problem fl2L fl3j . In particular, the interacting quintessence 
models of references 0, ^|, require the ratio of matter and dark energy densities to be constant at late times. 
The coupling between matter and quintessence is either motivated by high energy particle physics considerations Il2| 
or is constructed by requiring the final matter to dark energy ratio to be stable against perturbations [Til [Til Il6| . 
The nature of both DM and DE being unknown, there is no physical argument to exclude their interaction. On 
the contrary, arguments in favor of such interaction have been suggested |17|. and more recently they have been 
extended to include neutrinos 0. As a result of the interaction, the matter density drops with the scale factor 
a(t) of the Friedmann-Robertson-Walker metric more slowly than a~ 3 . The interacting quintessence model studied 
in the literature have been found to agree with observations of WMAP data and supernovae i but they require 
values of cosmological parameters different from those of WMAP (first year) concordance model. Observations of 
the large scale structure can also be used to constrain models. Recent data includes the SDSS 3 and 2dFGRS [j| 
measurements of the matter power spectrum. The analysis of 2dFGRS showed discrepancies with WMAP first year 
data but is in much closer agreement than SDSS with the results of WMAP third year data Q • 

Currently, there is no compelling evidence for DM-DE interaction [2fJ| and its (non-)existence must be decided 
observationally. It has been suggested that the skewness of the large scale matter distribution is a sensitive parameter 
to determine the difference in the clustering of baryons and dark matter resulting from the interaction [2l| . In this 
paper we shall study the effect of the interaction on the evolution of matter density perturbations during the radiation 
dominated period. The evolution of matter and radiation density perturbations provides powerful tools to constrain 
the physics of the dark sector j^. We will show how the shape of the matter power spectrum can be a directly 
related with the interaction and we will use the matter power spectrum measured by the 2-degree field galaxy redshift 
survey (2dFGRS, for short) to set constraints on the interaction. The outline of the paper is as follows: In Section 
II, we present a brief summary of the interacting quintessence model (IQM, hereafter). In Section III we describe the 
evolution of matter and radiation perturbations in models with dark matter and dark energy. In Section IV we discuss 
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some analytical solutions and in Section V we show how the slope of a scale invariant matter density perturbations 
has less power on small scales than non-interacting models. In Section VI we describe the results of Monte Carlo 
Markov chains constructed to compare the model with the observations. Finally, Section VII summarizes our main 
results and conclusions. 



II. THE INTERACTING QUINTESSENCE MODEL 



Most cosmological models implicitly assume that matter and dark energy interact only gravitationally. In the 
absence of an underlying symmetry that would suppress a matter - dark energy coupling (or interaction) there is no 
a priori reason for dismissing it. Cosmological models in which dark energy and matter do not evolve separately but 
interact with one another were first introduced to justify the small value of the cosmological constant |23|. Recently, 
various proposals at the fundamental level, including field Lagrangians, have been advanced to account for the coupling 
[2"i| . Scalar field Lagrangians coupled to matter generically do not generate scaling solutions with a long enough dark 
matter dominated period as required by structure formation . The phenomenological model we will be considering 
was constructed to account for late acceleration in the framework of Einstein relativity and to significantly alleviate 
the aforesaid coincidence problem IbH Ju& an d escapes the limits imposed by [25|. Here we shall describe its main 
features. For further details see Refs. {lltHiH^ . 

The model considers a spatially flat Friedmann-Robertson-Walker universe filled with radiation, baryons, dark 
matter (subscript, c) and dark energy (subscript, x). Its key assumption is that dark matter and dark energy are 
coupled by a term Q = 3Hc (p c + Px), that gauges the transfer of energy from the dark energy p x to the dark matter 
p c . The quantity c 2 is a small dimensionless constant that measures the strength of the interaction, and H — a -1 da/dt 
is the Hubble function. To satisfy the severe constraints imposed by local gravity experiments [ToL |27|. baryons and 
radiation couple to the other energy components only through gravity. Thus, the energy balance equations for dark 
matter and dark energy take the form 

^+3Hp c = Q, ^L+3H(l + w x )p x = -Q, (1) 
dt dt 

where w x — p x / p x < — 1/3 is the equation of state parameter of the dark energy fluid. 

Our ansatz for Q guarantees that the ratio between energy densities r = p c j ' p x tends to a fixed value at late times. 
This can be seen by studying the evolution of r which is governed by 

dr 2 (Pc + Px) 2 

— = -~3THr , T = -w x - c . (2) 

dt p c p x 

The stationary solutions of Eq. J2Jl follow from solving r s T(r s ) — 0. When w x is a constant these solutions are given 
by the roots of the quadratic expression 



rf = -l + 2b±2y/b(b-l), b = ~^>l. (3) 

As it can be checked by slightly perturbing Eq. J2J, the stationary rf solution is unstable while r~ is stable. The 
general solution of Eq. © can be written as 

r{x) = \ + + X ^ , (4) 

where x = (a/ao)~ A with A = 12 c 2 ^b(b — 1). In the range r~ < r < rf r(x) is a monotonic decreasing function. 
Thus, as the Universe expands, r(x) gently evolves from rf to the attractor solution r~ . The transition from one 
asymptotic solution to the other occurs only recently (see Fig. 2 of so we can take r ~ rf during a fairly large 
part of the history of the Universe. Finally, the constraint on rf ~ const implies that r and c 2 are not independent, 
but linked by c 2 (rf + l) 2 = rf | w x |, so the product c 2 rf ~| w x is of order unity. 

We would like to remark that the above ansatz for Q is not arbitrary. It was chosen so that the ratio between dark 
matter and dark energy densities tends to a fixed value at late times, thereby alleviating the coincidence problem 
[ill Il5l |26| . It also yields a constant but unstable ratio at early times. It is hard to imagine a simpler expression 
for Q entailing these two key properties. Likewise, it is only fair to acknowledge that the aforesaid expression can be 
re-interpreted as implying, at late times, an effective exponential potential for the quintessence field. This well-known 
result was derived by Zimdahl et al. Likewise, in we remarked that the effective potential of the IQM model 
exhibits a power-law dependence on the quintessence field at early times and an exponential dependence at late times. 
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Near r « r s the balance Eqs. (JIJ can be approximated by 

1 d Pc l-c 2 (l + l/r7) 



p c dr ' c 2 (rf ~ rj)(r - rj) 
J_ dfx_ ^ 1 + Wx + c 2 (l + r£ 
p K dr ' c 2 (rf - rs){r - n) 

For w x ~ constant, these equations can be integrated to 



(5) 



Peocfl- 3 ! 1 - 2 * 1 ^)], ^oca-^+^+^+'^L (6) 

Notice that the condition r(r\r) = implies that the exponents in the energy densities, Eq. (JHJ), coincide. Interestingly, 
these results are not only valid when the dark energy is a quintessence field (i.e., —1 < w x < —1/3), they also apply 
when the dark energy is of phantom type (i.e., w x < —1), either a scalar field with the "wrong sign" for the kinetic 
energy term, a fe-essence field, or a tachyon field jllj . 

Near the attractor, dark matter and dark energy dominate the expansion and Fricdmann equation becomes simply 
3 H 2 — n(p c + p x ) and a oc f( 2 / 3 )l 1 + w ^+ c ( 1 + r s )] . The results presented here significantly alleviates the coincidence 
problem but they do not solve it in full. For this purpose, one needs to show that the attractor was reached only 
recently -or that we are very close to it- and that r~ is of order unity. In fact, the value of r~ cannot be derived from 
data and must be understood as an input parameter. This is also the case of a handful of key cosmic quantities such 
as the current value of the cosmic background temperature, the Hubble constant or the ratio between the number of 
baryons and photons. 



III. LINEAR PERTURBATIONS 



As the scalar field is coupled just to dark matter and since dark matter and quintessence are coupled to baryons and 
photons only gravitationally, there is no transfer of energy or momentum from the scalar field to baryons or radiation 
and their evolution is the same as in non-interacting models. In the synchronous gauge and for a flat space-time, the 
line element is given by: ds 2 — a 2 (r)[— dr 2 + (5k, ij + hij)dx l dx^] 1 where r is the conformal time, a the scale factor 
and &K,ii is Kronecker's delta tensor. Only two functions h and p are necessary to characterize the scalar mode of 
the metric perturbations hij |28j . Assuming the dark energy energy-momentum tensor is free of anisotropic stresses, 
the equations describing the dark matter and dark energy evolution in the synchronous gauge are: 

4 = -(l + w x )(9 x + ^)-3H(l-w x )5 x -m 2 (l-w 2 x ) ( ^+ 3Hc 2 {d x +rd c ) 1 (7) 

9 X = 2W x + -^5 x -3H-^—(l + r)9 x , (8) 
1 + w x l + w x 

h S 

S c = -e c -^-3Hc 2 {5 c + —) : (9) 
9 C = -H9 c + 3Hc 2 {l + l/r)9 x . (10) 

In these expressions, the derivatives are taken with respect to the conformal time t, S is the density fluctuation, 9 
the divergence of the fluid velocity, h the gravitational potential, k the wavenumber of a Fourier mode, TL = a/ a and, 
r is the ratio of the background cold dark matter to the dark energy density. We also assume that the dark energy 
has constant equation of state parameter w x = constant and sound speed c s , x = 1. In this gauge, (SP/Sp) x 5 x — 
c 2 s x 8 x + 3H9 X (1 + w x )(c 2 s x - w x )/k 2 . As noted in [13, Eq. (JTUJl was mistyped in (l9|. 

Equations l(7l)- (|10|) do not form a closed set. They must be supplemented with the equations describing the evolution 
of the coupled baryon-photon fluid, neutrinos and gravitational potentials. For the potentials, the only relevant 
quantity is the trace of the metric perturbation, h. Its time evolution can be derived from Einstein's equations: 

h + Hk = -m 2 (i + 3 (sp/6 P ) 2 s l ) s.n, , (11) 

i 

where the sum is over all matter fluids and scalar fields; fi^ is the energy density of fluid i in units of the critical 
density. With respect to baryons, photons and neutrinos, they interact with the DM and DE only through gravity. 
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Figure 1: Evolution of the gravitational potential (upper panels) and the cold dark matter density perturbations (lower panels) 
for three modes: k — 0.01 (left) k = 0.1 (center) and k — l/i^Mpc (right panels). We study the evolution of each mode in 
three different cosmological models: the concordance model (c 2 = 0, solid line) and two interacting quintessence models with 
the same cosmological parameters, c 2 = 10 -3 (dotted) and c 2 = 6 x 10 -3 (dashed line). 



The coupled evolution of dark matter, dark energy, baryon, photon and, optionally, neutrino density perturbations 
and gravitational fields can not be solved analytically. To compute numerically the solution, we have implemented 
Eqs. (f7f> — (fTTf) into the CMBFAST code [2!j ■ In Fig. Q]we show the evolution of the potential, h, and the matter 
density perturbation, S c , for three modes of wavelength k = 0.01,0.1 and lh Mpc -1 and for three different values of 
the DE decay rate: c 2 = (solid line), c 2 = 10~ 3 (dotted line) and c 2 = 6 x 10~ 3 (dashed line). In all the cases, the 
cosmological parameters defining the background model are: f2 c = 0.26, = 0.04, Q, x = 0.7, w x = —0.9 and the 
Hubble constant Ho = 70 km/s Mpc -1 . 

As panels ^3 and ^ illustrate, modes that enter the horizon before matter-radiation equality grow slower with 
increasing interaction rate. As a result, the matter power spectrum on those scales will have less power than in non- 
interacting models. To obtain some insight on the behavior on the evolution of matter density perturbations, we will 
be considering some limiting cases where analytic solutions exist. For simplicity, we shall assume the dynamical effect 
of baryons and neutrinos in the evolution of dark matter and dark energy density perturbations can be neglected. 
The result of combining Eqs. © and IjlOfl. considering only the leading terms at first order in the interaction c 2 , is 



I + H(l + 3c 2 )S c + 3c*H 2 + ^P^r ) S c = -\(h + Hh) + c 2 F(S Xl S x ), 

\ 2 (l + w x )AJ 2 



(12) 



with 



F(S X ,S X ) 
A 



3H g x + 3 2fi \Hr + l)(l-w x -c 2 ) l-3w B 



(1 + w x )A 



i + (m/k) 2 {i-w 2 x ). 



(13) 
(14) 



Equation l|12|) corresponds to a damped harmonic oscillator (with real or imaginary frequencies) with a forcing term 
that, since c 2 <C 1, is dominated by the time evolution of the gravitational potential. In this approximation, Eq. i|ll|) 
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gives 



Q c f 3(l-w x )(3H/k) 2 ^\ c2 {3c 2 (r + l) t l-3w B \ 



2{ i + A C )- C \SVT^)A + —^) 



In the limit c 2 = this equation coincides with the evolution of matter density perturbation in non-interacting 
cosmologies. The effect of the interaction is to increase the friction term and to modify the oscillation frequency. 
The term in square brackets accounts for the self attractive force acting on the perturbation; the extra contribution 
arises due to the interaction. The effect of the DM-DE coupling can be understood as a modification of the effective 
gravitational constant. This result was previously found in |30| and |3l| . both models with a different interaction 
ansatz. In our case, the interaction provides a new physical effect not present in other models: if w x ~ —1, the 
second term in the square parenthesis could dominate and matter density perturbations would stop growing and start 
oscillating. 

IV. EVOLUTION OF MATTER DENSITY PERTURBATIONS 
A. Superhorizon sized perturbations 

The time variation of dark matter and dark energy densities have analytic expressions in terms of the expansion 
factor (Eq. ©). In terms of the time variable log a, analytic solutions can be found for the evolution of superhorizon 
sized perturbations. Using this new time variable Eq. (|15(l can be written as 

K + (^f^ + 3c 2 ) S' c + 3c 2 (i^) S c = \(1 + 3< s )^ , (16) 

where prime denotes derivatives with respect to log a. The subindex B stands for background quantities. The behavior 
of Sps/pB at scales larger than the Jean's length can be parameterized as: 



Spb _ f 8pb \ f _a_ N '' " 



Pb \ Pb J Hi \ a-H! 



(17) 



where {5pbPb)h, is the amplitude of the mode under consideration at horizon crossing, at time an r - After a brief 
transient period, the evolution of the dark matter and dark energy perturbations will be given by the inhomogencous 
solution associated with the time evolution of the gravitational potential: S c ~ h ~ a p / 2 . The solutions of eq (|16fl are: 
P = (4; 2 — 6.6c 2 ) in the radiation and cold dark matter dominated periods, respectively. For non-interacting models, 
the well known solutions are p — (4; 2). These solutions were to be expected; as discussed in |32|], if p c ~ fl ~°S 
being a a constant, then p ~ 2{a — 2). In the radiation epoch, a = 4, and in the matter epoch a = 3(1 — c 2 ). 
Thus, during matter domination, the growth of dark matter density perturbations slows down with respect to those 
of non- interacting models but, in general, the evolution of superhorizon sized perturbations is not significantly altered 
by the interaction. 

B. Subhorizon sized matter perturbations 

For perturbations inside the horizon, Eqs. Q and JHJ) have the approximate solution, 

Sx = 3 {j) c2rSc - (18) 

Since c 2 <c 1, we have that fl c S c 3> £l x $x and the force term in Eq. (|15f) is dominated by the perturbations in the 
photons field during the radiation epoch. In the small scale limit (37i/fc <C 1), A ~ 1 and 



H(l + 3c 2 )<5 c - AirGa 2 p c 



3c 2 rw x 



ATtG Pc a?{r + 1)(1 + w x ) 



H 2 rJ 7 (S 7 . (19) 



Even this simplified equation does not have simple analytic solutions. The slower growth of matter density pertur- 
bations in the IQM compared with non-interacting ones can be understood analyzing the different coefficients: (A) 
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the interaction increases the friction term, damping more rapidly the homogeneous solution and (B) it decreases the 
gravitational force acting on the perturbation. At very early times, when fi c < 6c 2 | w x /{\ + w x ) , and well within 
the horizon perturbations on the photon field oscillate and matter perturbations do not grow but undergo damped 
oscillations [33. The characteristic time-scale of the growth of matter density perturbations is the mean free-fall 
time, tjf ~ (Gp c ) -1 / 2 . During the radiation dominated regime the expansion rate is fixed by Friedmann equation: 
H ~ */ Gp 1 and since tff >> H^ 1 , perturbations only grow logarithmically, not as a power law. In our IQM this 
effect is more severe. First, at all times the matter density is smaller than in non-interacting models but with the same 
values for the cosmological parameters today. Second, the effective gravitational force is reduced. Thus, the mean 
free-fall time increases and density perturbations grow slower (or get even erased), compared with a non-interacting 
model. 

Scalar fields coupled to matter would modify gravity inducing an extra attractive force. A repulsive effect could be 
obtained by the exchange of vector bosons. It was first suggested |34| that phantom scalar fields with non-standard 
kinetic term coupled to matter would give rise to a long-range repulsive force. In our phenomenological model, the 
decrease of the gravitational coupling in Eq. I|19(l is due to our specific ansatz Q for the dark matter - dark energy 
interaction. 



C. Comparison with other interacting models 

Interacting quintessence models couple dark matter and dark energy so the en ergy momentum tensor of the DM 
and DE are not separately conserved but obey {T^ v + T^ L )- I _ L = 0. In references 0,123 the coupling is chosen such 

that p c + 3Hp c — (167TG/3) 1 / 2 (3p c 4> x , where (f> x is the scalar field describing the dark energy component and (3 the 
decay rate coefficient whom, in general, is a time varying function. By assuming that the scalar field couples to dark 
matter only, the evolution of matter density perturbations in the synchronous gauge is given by 

5c = -0 c -±h-±{fi V ), (20) 

6 C = -H(l-2(3x)0 c -2f3k 2 (j). (21) 

In this expression, ip = (4-ttG /2>) l / 2 5<f> x is the perturbation in the scalar field and x its kinetic energy. The evolution 
of the gravitational field does not depend on the specific interaction ansatz, and is given by Eq. iJUJ. For subhorizon 
sized perturbations, ip — fc~ 2 7i(/3f2 c 5 c ) and, in the radiation dominated regime, matter perturbations evolve as 

5 C + H[l - 2/3x]S c - 4vrGa 2 [l + \p 2 ]p c 5 c = ^?i 2 2fi 7 5 7 . (22) 

As discussed above, during the radiation period the background expansion rate is fixed by Friedmann's equation, but 
the mean free-fall time is now: tff ~ [G(l + 4/3 2 /3)p c ] 1 / 2 . Due to the interaction, the dark matter density at any given 
time is smaller than in a non-interacting model with the same cosmological parameters, and the difference increases 
with p. Likewise, the effective gravitational constant increases, but since the dependence is second order in /3, one 
would expect tff to be smaller than in non-interacting models. This statement depends on the particular interaction 
ansatz. Since perturbations evolve as if the Newton's gravitational constant was a factor (1 + 4/3 2 /3) larger, the 
interaction with the scalar field could make density perturbations to grow faster during the matter dominated regime 
due to a larger local gravity. This effect could compensate the slow growth during the radiation dominated regime 
and enhancing the clustering of dark matter perturbations compared with the uncoupled case, as found in |3l]|. But 
even in this case, the amplitude of the matter power spectrum was smaller in the range (0.01 — 0A)h Mpc" 1 . 

V. THE EFFECT OF THE INTERACTION ON THE MATTER POWER SPECTRUM 

In the previous section we have shown that the interaction slows the growth of matter density perturbations. Only 
the slower growth of perturbations in the radiation dominated regime will have a significant impact on the matter power 
spectrum today. For comparison, we shall assume that in interacting and non-interacting models density perturbations 
have the same amplitude when they come within the horizon. For non-interacting models, this prescription leads to 
the so called Harrison-Zeldovich power spectrum [35|, characterized by a functional form P(k) ~ fc" s with n s = 1 
on large scales. During the matter epoch, if density perturbations evolve with the scale factor as S c ~ dpi 2 and the 
background energy density as p c ~ aT a (with a =const) the power spectrum will scale with wavenumber as 



P(k) ~ fc- 3+2 PA"- 2 ) 



(23) 
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Figure 2: (a) Matter power spectra for the interacting quintessence model with different rates of energy transfer. From top to 
bottom c 2 = 0, 10~ 3 , 6 x 10~ 3 , 10~ 2 . We took the present value of cosmological parameters to be: fib = 0.04, fi c dm = 0.23, 
Ho = 72km s~ 1 /Mpc, Qa = 0.73 the dark energy equation of state w x — —0.9 and the slope of the matter power spectrum 
at large scales n 3 = 1. (b) The same for mixed dark matter models with one single species of massive neutrinos. From top to 
bottom, the fraction of energy density in the form of neutrinos is: Q. v = 0.01,0.05,0.1,0.2. As before, the total dark matter 
energy density was Q,dm = 0.23, the rest was cold dark matter, (c) Variation of the slope of P(k) with c 2 , and (d) with massive 
neutrinos. The slope was computed from a straight line fit to the data in the interval k — [0.1, l]ft _1 Mpc. 



During the matter dominated period p ~ 2(ce — 2) + 0.6c 2 and the slope of a scale-invariant spectrum is n s — 1, with 
a very weak dependence on the interaction. 

The slope of the matter power spectrum on scales k > k eq is determined by the growth rate of subhorizon sized 
matter perturbations during radiation domination. If a mode that crosses the horizon before matter-radiation equality 
(k > k eq ) grows as S c ~ T q ' 2 during the radiation dominated period, then the amplitude of the power spectrum today 
would be P(k) — P(k eq )(k eq / k)~ 3+q . For cold dark matter models, dark matter perturbations experience only a 
logarithmic growth, so models with less growth will have less power at small scales as do, for example, mixed dark 
matter models 36], i.e., models containing a significant fraction of massive neutrinos. 

In Fig. we plot the power spectrum for different interacting quintessence models. All models have the cosmo- 
logical parameters of the WMAP first year concordance model |2j- From top to bottom, c 2 = 0, 10~ 4 , 1CP 3 , 10~ 2 ; the 
normalization is arbitrary. Similarly, in Fig. ^jp we plot the matter power spectrum of mixed dark matter models with 
one species of massive neutrinos, for different neutrino masses: m u — 0,0.1, 1, 10 eV. With increasing decay rate or 
neutrino mass, the matter power spectrum shows larger oscillations, due to the increased ratio of baryons to dark mat- 
ter. The slope decreases with increasing c 2 and m„. As explained above, potential wells are shallower with increasing 
c 2 ; matter perturbations during radiation domination are damped similarly as do in models with massive neutrinos. 
In Fig. |2t we plot the change in the slope of the matter power spectrum as a function of the energy transfer rate and 
in|5p as a function of the neutrino mass. As the slope changes smoothly from large to small scales, for convenience we 
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computed the slope as a straight line fit to the data in the interval k = [0.1, 1.0]/i _1 Mpc. In both cases the behavior 
is rather similar: for low values of neutrino mass and interaction coupling, the slope is approximately —2 and roughly 
constant. When parameters are increased in either model, the slope decreases. Observations of large scale structure 
that constrain the neutrino mass can also be used to set constraints on the strength of DM-DE coupling during 
the radiation dominated period. These constraints are complementary to those coming from skewness of the matter 
density field, that are sensitive to the interaction at much lower redshifts |2l|. Fig. [2] shows a significant difference 
between massive neutrinos and interacting quintessence: in IQM the maximum of the matter power spectrum shifts 
to larger scales. At larger c 2 , the dark matter density becomes smaller at any given redshift and the matter radiation 
equality is delayed. This does not happen in models with massive neutrinos where matter-radiation equality occurs 
always at the same redshift. 



VI. OBSERVATIONAL CONSTRAINTS ON DARK MATTER - DARK ENERGY COUPLING 

Since the interaction affects the slope of the matter power spectrum, we used the 2dFGRS data [j| to constrain 
c 2 . We used a Monte Carlo Markov chain to run the CMBFAST code, adapted to solve the IQM described above, 
through a 7-dimensional parameter space: (A, fl(,h 2 , Q c h 2 , Hq, n Sl c 2 , w x ) where A is the normalization of the matter 
power spectrum, f2&, Q c are the baryon and cold dark matter fraction in units of the critical density, n s is the slope 
of the matter power spectrum at large scales, c 2 measures the transfer rate of dark energy into dark matter and w x 
is the dark energy equation of state parameter. Hereafter Hq = 100ft km/s Mpc -1 . It is common practice to call h 
the Hubble constant in units of 100km s _1 /Mpc and we shall follow this convention. It should not be confused with 
the gravitational potential in the synchronous gauge introduced in Sec. III. We limit our study to flat models, so the 
fraction of dark energy is fixed by the Friedmann equation il x + fit, + il c + f2 7 = 1, where f2 7 is the photon energy 
density, and all densities are measured in units of the critical density. To simplify, we studied only adiabatic initial 
conditions and initial power spectrum with no running on the spectral index. We did not include reionization, or 
gravitational waves, since they have little effect on the matter power spectrum. 

Since we are interested in constraining c 2 from the shape of the matter power spectrum, we have to correct for 
non-linear effects. We followed the 2dFGRS team and assumed the non-linear biasing to be well described by 

Pgal{k) = 6 2 ^ fc * fWm(*0- (24) 

We used A g = 1.4 and Q = 4.6. We marginalized over the bias factor b. We did not use the SDSS galaxy power 
spectrum because these data were analyzed in real space where non-linear effects are more important. We used 
the likelihood codes provided by the 2dFGRS team [5|. As priors, we imposed our chains to take values within 
the intervals: A = [0.5,2.0] in units of COBE normalization, h = [0.4,1.1], tt b h 2 = [0.00,0.05], tt c h 2 = [0.0,0.5], 
n s — [0.80, 1.2], w x = [—0.5, —1.0] and c 2 = [0,0.05]. We run the chain for 10 5 models, that were sufficient to reach 
convergence. In Fig. [3]we plot the marginalized likelihood function obtained from the posterior distribution of models. 
The likelihood is very non-gaussian, reflecting the fact that models do not depend linearly on c 2 . The data are rather 
insensitive to c 2 < 10~ 3 since the slope does not change significantly up to that value (see Fig. As discussed 

above, increasing the interaction rate leads to smaller fraction of dark matter during the radiation dominated period 
and shallower potential wells, larger free-fall times and, as a result, the amplitude of the matter power spectrum is 
damped (see Fig. QJ). 

In Fig. 0]we show the joint confidence contours at the 68%, 95% and 99.99% level for pairs of parameters after 
marginalizing over the rest. Our prior on the spectral index was too restrictive and did not allow the chains to sample 
all the parameter space allowed by the data. Therefore, we can not draw definitive conclusions about the confidence 
intervals for all of the parameters. The figure does show that the data at present do not have enough statistical power 
to discriminate the IQM from non-interacting ones. Our la confidence levels and upper limits for the cosmological 
parameters are c 2 < 3 x 10 -3 , Q c h 2 = 0.1 ± 0.02, H = 83l^ km s _1 /Mpc. The data are rather insensitive to w x 
and baryon fraction. Models with c 2 = are compatible with the 2dFGRS data at the la level. The data show a full 
degeneracy with respect to c 2 up to c ~ 10 -3 , in contrast with the results of [19( obtained using the 1st year WMAP 
data. There the data preferred interacting quintessence models with respect to non-interacting ones but this was an 
artifact of our parameter space since we restricted the normalization to be that of COBE, penalizing the concordance 
model that prefers a lower normalization. A full discussion including WMAP 3rd year data will be deferred to a 
forthcoming paper. 
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Figure 3: Marginalized likelihood function for the 2dFGRS data. 
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Figure 4: Marginalized likelihood function for the 2dFGRS data. 
VII. DISCUSSION 

Interacting quintessence models have been constructed to solve the coincidence problem. They make specific 
predictions that can be checked against observations of large scale structure. In this paper we have shown that the 
interaction induces measurable effects in the growth of the matter density perturbations and modifies the power 
spectrum on small scales. To summarize, if dark energy decays into dark matter the background model has less dark 
matter and more dark energy in the past compared to non-decaying models. Since the dark energy does not cluster 
on small scales it does not contribute to the growth of density perturbations. The mean free-fall time increases and 
perturbations grow slower than in non-interacting models. The slower growth of matter density perturbations during 
the radiation and matter period, results on a damping of the matter power spectrum on those scales that cross the 
horizon before matter-radiation equality, but it does not change the slope on large scales. 

The combined effect of shifting the scale of matter radiation equality and changing the slope of matter power 
spectrum at small scales is a distinctive feature of interacting models where the dark energy does not cluster on 
small scales. Measurements of matter power spectrum could eventually reach enough statistical power to discriminate 
between interactin g an d non-interacting models. For example, the spectrum obtained from Lymann-a absorption lines 
on quasar spectra 37] probe the matter power spectrum at redshifts in the interval (2,4), where non- linear evolution 
has not yet erased the primordial information down to megaparsec scale. The use of more precise information on 
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small and large scales, could set tighter bounds on the interaction of dark matter and dark energy. 

Our previous results 0] and those presented here indicate that the IQM fits the observational data as well as non- 
interacting models, alleviates the coincidence problem and provides a unified picture of dark matter and dark energy. 
It predicts a damping on the matter power spectrum on small scales that can be used, together with the delay on 
the matter-radiation equality, to discriminate it from non-interacting models. The slower growth of subhorizon sized 
matter density perturbations within the horizon provides a clean observational test to proof or rule out a DM-DE 
coupling. 
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